The stability and the shape of the heaviest nuclei 
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In this paper, we report a systematic study of the heaviest nuclei within the relativistic mean 
field (RMF) model. By comparing our results with those of the Hartree-Fock-Bogoliubov method 
(HFB) and the finite range droplet model (FRDM), the stability and the shape of the heaviest 
nuclei are discussed. The theoretical predictions as well as the existing experimental data indicate 
that the experimentally synthesized superheavy nuclei are in between the fission stability line, the 
line connecting the nucleus with maximum binding energy per nucleon in each isotopic chain, and 
the /3-stability line, the line connecting the nucleus with maximum binding energy per nucleon in 
each isobaric chain. It is shown that both the fission stability line and the /3-stability line tend to 
be more proton rich in the superheavy region. Meanwhile, all the three theoretical models predict 
most synthesized superheavy nuclei to be deformed. 

PACS numbers: 21.10.Dr, 27.90.+b, 21.60.-n 



I. INTRODUCTION 

Recent claims of successful syntheses of superheavy el- 
ements 115 and 113 0,0 have aroused new enthusiasm 
about studies of superheavy nuclei in the nuclear physics 
community (see Refs. 0, 0, and references therein). 
For a review of recent experimental progress of this sub- 
ject, we refer the reader to Refs. [a, Q. Conventional 
liquid drop models of finite nuclei forbid the existence 
of any nuclei with a proton number larger than 100, i.e. 
superheavy nuclei, due to the destructive Coulomb force. 
However, shell effects are found to be able to stabilize 
these nuclei, and therefore explain their very existences 
0, 0. It has long been predicted that there exist a 
large number of relatively long-lived superheavy nuclei, 
the so-called superheavy island, which is separated in 
neutron and proton numbers from the known heavy ele- 
ments by a region of much higher instability. Although 
the experimentally-synthesized superheavy nuclei are in- 
deed very heavy, it is generally believed that they are not 
examples of the originally sought island of superheavy el- 
ements. 

On the theoretical side, a lot of efforts have been made 
to interpret the experimental results and make various 
predictions. A short review of the theoretical activities 
can be found in Ref. Nowadays, there are several cat- 
egories of theoretical models often used to study super- 
heavy nuclei: The first category is the liquid drop model 
and its many variants, such as the finite-range droplet 
model (FRDM) [13; The second category is the non- 
relativistic Skyme-Hartree-Fock model; The third cate- 
gory is the relativistic mean field model. Using these 
models, fission barriers, alpha-decay energies, shell clo- 
sures, single-particle spectra, and so on have been exten- 
sively discussed pj. 

In recent years, the relativistic mean field model has 
received much attention due to its natural description 
of the spin-orbit interaction, the saturation properties of 



symmetric nuclear matter, and many other things that 
non-relativistic models have some difficulties to explain 
|llj. It has also been widely employed to study super- 
heavy nuclei (a short review can be found in Ref. 5]). 
Due to the amount of computer resources needed, studies 
of superheavy nuclei in the relativistic mean field model 
have often been limited to either spherical assumption 
or a small part of the superheavy region. In this pa- 
per, we report the first systematic study of superheavy 
nuclei within the relativistic mean field model with the 
deformation effect and the pairing correlation properly 
treated. By comparing our RMF+BCS calculations with 
those of the Hartree-Fock-Bogoliubov (HFB) and finite 
range droplet model (FRDM), we hope to obtain some 
hints for the search of the superheavy island. 

In particular, we would like to adress two interesting 
subjects: the stability and the shape of the heaviest nu- 
clei. The stability of a superheavy nucleus is a very sub- 
tle subject. It is determined by many competing decay 
modes: alpha decay, spontaneous fission, beta decay, and 
etc. Except beta decay, the other two are very difficult 
to describe quantitively, which would involve complicated 
lifetime calculations. Beta decay, on the other hand, can 
be understood much more easily from the energy point 
of view, i.e. nuclei near the /3-stability line are stable 
against beta decay. Similarly, a useful concept is the fis- 
sion stability line, the line connecting the nucleus with 
maximum binding energy per nucleon in each isotopic 
chain, which is related to the minimum Q value of fis- 
sion [12 ]. At this line, nuclei with fixed proton num- 
ber have maximum binding energy per nucleon; therefore 
they would be stable against neutron emission, which can 
play an important role to synthesize superheavy nuclei at 
the first place 0,0. Hence, one would expect superheavy 
nuclei should not deviate too far from this line. 

The shape of the heaviest nucleus can also influence its 
stability greatly 01 ■ The original island of superheavy 
elements is predicted to be around Z = 114 and N = 184 
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[Tij mainly due to the fact that 184114 is predicted to be 
a doubly magic nucleus with spherical shape, where shell 
effect is the strongest. Recent investigations have pro- 
vided somehow conflicting predictions for the next dou- 
bly magic system, for example, the non-relativistic forces 
SkM* and SkP seem to prefer 184126 instead of 184114 
|l5j . A more complete summary of various predictions 
of different models can be found in Ref. [5j- In this 
paper, we would like to compare our RMF+BCS predic- 
tions with those of the HFB model and the FRDM model 
in order to see whether a large number of spherical nuclei 
exist in these models, which would indicate the existence 
or nonexistence of the next doubly magic nucleus, or (less 
ambitiously) the next neutron (proton) shell closure. 

This paper is organized as follows. In Section II, we 
briefly introduce the relativistic mean field model and ex- 
plain the numerical details of our calculation. In Section 
III, the stabilities of superheavy nuclei are studied from 
the energy point of view. In Section IV, the shapes of su- 
perheavy nuclei predicted by different theoretical models 
are compared. The whole work is summarized in Section 
V. 



II. THEORETICAL FRAMEWORK 

In this section, we briefly describe the RMF+BCS cal- 
culations. The RMF+BCS calculations have been car- 
ried out using the model Lagrangian density with non- 
linear terms for both a and uj mesons as described in 
detail in Refs. 0,0, which is given by 

C = ^(n^ - M)if) 

+ \d^ad^a - \m 2 a a 2 - ^92<J 3 - ^93^ - g^wp 

- tn^VL^ + \mlu)^ + \g^^f - g^jfipcj/, 

- \R% v R ailv + \m 2 p pl P ^ - g^T^p^ 

(1) 

where all symbols have their usual meanings. The corre- 
sponding Dirac equation for nucleons and Klein-Gordon 
equations for mesons obtained with the mean-field ap- 
proximation and the no-sea approximation are solved by 
the expansion method on the axially deformed harmonic- 
oscillator basis [Isj . The number of shells used for ex- 
panding the nucleon and meson wave functions is cho- 
sen as Nf = Nb = 20. More shells have been tested 
for convergence considerations. Quadrupole constrained 
calculations 0] have been performed for all the nuclei 
considered here in order to obtain their energy surfaces 
and determine the corresponding ground-state deforma- 
tions. 

The pairing correlation plays an important role in stud- 
ies of open-shell nuclei. It is also true for superheavy nu- 
clei. In Ref. , it was shown that the use of a zero-range 
(5-force in the particle-particle channel can bring down 
the fission barrier compared to the use of a constant-gap 
pairing method. In the present calculation, the pairing 



correlation is treated by a state-dependent BCS method 
ppf. More specifically, the pairing force used is of the 
volume type 

V = V 6(r-7). (2) 

In the past years, whether the pairing correlation in fi- 
nite nuclei is of a volume type or surface type has been 
discussed a lot , but it seems that more investigations are 
still needed to reach a definite conclusion [2l|; therefore, 
to limit the number of free parameters, we do not in- 
troduce explicitly any density dependence. On the other 
hand, to describe simultaneously both light and heavy 
nuclei |l6l |. we introduce a weak mass number depen- 
dence to the pairing strength, i.e. 

V a = 300+ 120/ A 1/3 , (3) 

which is purely phenomenological except for the A 1 / 3 de- 
pendence |22j |. For nuclei with an odd- number of nucle- 
ons, the blocking effect has been treated within the BCS 
framework 0,0. A more detailed description of the pair- 
ing method can be found in Ref. Q. 

In the mean field channel, the effective force TMA 17] 
is used. The effective force TMA was first proposed to 
describe simultaneously both light and heavy nuclei. It 
in fact originated from two other very successful param- 
eter sets: TM1 and TM2. TM1 aimed to describe the 
ground-state properties of heavy nuclei (^4 > 40) and 
TM2 those of light nuclei (A < 40). On the one hand, 
TMA inherited TMl's favorable property of being able 
to reproduce the essential feature of the equation of state 
and the vector and the scalar self-energies of the relativis- 
tic Bruckner-Hartree-Fock theory for nuclear matter [2j| . 
On the other hand, it can also describe light nuclei very 
well. Its success comes from a weak mass dependence, 
which smoothly interpolates the TM1 and the TM2 pa- 
rameter sets. We may interpret this mass dependence 
as a mean to effectively express the quantum fluctua- 
tions beyond the mean field level and/or the softness of 
the nuclear ground states in deformation, pairing and 
alpha clustering in light nuclei. Compared to other suc- 
cessful effective forces, such as NL3, the description of 
finite nuclei is of similar quality or slightly better (lt| . 
but TMA(TMl) yields a much softer equation of state 
at high density, which seems to be favored by current 
experimental results. 

Finally, we would like to mention the strategy we used 
to confine our study to a reasonable number of nuclei 
since calculations of superheavy nuclei cost a lot of time. 
The number of nuclei is determined by including all of 
those compiled in Ref. [24( and extending the corre- 
sponding proton-rich limit and neutron-rich limit of a 
certain isotopic chain by ten more nuclei, respec- 
tively (see also Fig|2J) . Using such a strategy, the number 
of superheavy nuclei investigated in the present study is 
around 600. 
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FIG. 1: (color online) a-decay energies of the 278 113 a-decay 
chain. The latest experimental results of Morita et al. @] are 
compared with the predictions of the RMF+BCS calculation 
and those compiled in Audi03 (obtained using "systematic 
trends") 0. 



III. THE STABILITY OF SUPERHEAVY 
ELEMENTS 



Before we study the stability of the heaviest nuclei, 
it is necessary to stress that although the parameters 
of the RMF+BCS model have not been constrained by 
any information from the superheavy region, the agree- 
ment of its predictions with existing experimental data 
is quite good in general, which has been demonstrated in 
the entire region (see Refs. 0, 0, |2^, |2^, |27j and refer- 
ences therein). Recently, Morita et al. has reported the 
synthesis of the new element 113 2]. In Fig. ^ the ex- 
perimental a-decay energies of the 278 113 a-decay chain 
are compared with our predictions and those compiled 
in Audi03 (obtained using "systematic trends" [2J] . It is 
clearly seen that the agreement is remarkable. A more 
detailed study is underway and will be reported some- 
where else. 

It was argued from the energy point of view that the 
fission stability line, the line connecting the nucleus with 
maximum binding energy per nucleon in each isotopic 
chain, plays an important role in studies of the stabil- 
ity of the heaviest nucleus in Ref. ^2] • I n Fig. |21 the 
binding energies per nucleon of the 600 nuclei we calcu- 
lated are plotted as functions of the neutron number iV 
(upper panel). The experimental data (lower panel) are 
taken from Ref. |24|. The nucleus with maximum bind- 
ing energy per nucleon in each isotopic chain is denoted 
by a triangle. Those nuclei in different isotopic chains 
are then connected by dashed lines. It can be clearly 
seen that the area we investigated has in fact included 
the most bound nucleus in each isotopic chain. The the- 
oretical curves are almost identical to the experimental 
ones (see also Fig. |3J). Therefore, just as expected, all 
the experimental syntheses are indeed near the fission 
stability line [ll. Since the results of the HFB-8 [H 
and FRDM |l0j mass formulae are very similar to our 
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FIG. 2: (color online). Binding energies per nucleon of super- 
heavy nuclei with Z — 101-120 and N = 134-190 as functions 
of the neutron number N. The theoretical predictions are 
compared with existing experimental data (including those 
obtained using "systematic trends") |24|. Different isotopes 
are connected by solid lines and ordered from top to bot- 
tom with increasing Z. The nucleus with maximum binding 
energy per nucleon in each isotopic chain is denoted by a tri- 
angle, and those of different isotopic chains are connected by 
dashed lines. 



calculations, they are not shown in this figure. Here, a 
few words about HFB-8 and FRDM are in place. The 
mass table of the finite-range droplet model has been 
around for more than ten years 10]. By carefully ad- 
justing its parameters (about thirty) to the saturation 
properties of nuclear matter and the binding energies 
of around 1000 nuclei, it obtained a root-mean-square 
deviation of about 0.6 MeV for the binding energies of 
all the experimentally-known nuclei. HFB-8 is a rather 
new mass table based on the Hartree-Fock-Bogoliubov 
method |28j |. By performing particle number projection, 
incorporating phenomenologically both the Wigner en- 
ergy and the rotational energy, and meanwhile adjusting 
its parameters (around 20) to fit the saturation properties 
of nuclear matter and the binding energies of about 2000 
nuclei, it achieved a similar quality to that of the FRDM 
model. More comparisons between these two models and 
current RMF models for nuclear ground-state properties 
can be found in Ref. 1291. 
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FIG. 3: (color online). The fission stability line and the /3-stability line as functions of neutron number N and proton number Z. 
The theoretical predictions (RMF+BCS, HFB-8 |U and FRDM [HI) are compared with existing experimental data (inclu ding 
those obtained using "systematic trends") [2^1 . The smooth dash-dot-dot line is the /3-stability line of Bohr-Mottelson |13| 
(see text). The small black dots represent experimentally synthesized superheavy nuclei (SHE) (including those obtained using 
"systematic trends") |24| . 



In Fig. |21 the fission stability line and the /3-stability 
line are plotted as functions of N and Z. The three the- 
oretical predictions (RMF+BCS, HFB-8 and FRDM 
[To|) are compared with existing experimental data pjj . 
For reference, the phcnomcnological /3-stability line 

N - Z = 6.0 x 10- 3 A 5/3 (4) 

derived from the Bethe-Weizsacker mass formula in Ref. 
[R| is also shown. Several interesting things can be 
learned immediately from Fig. [21 First, the microscopic 
(RMF+BCS, HFB-8 and FRDM) /3-stability lines agree 
with each other very well, but they are bent up a lit- 
tle bit compared to the phenomenological /3-stability line 
of Bohr-Mottelson. This difference can slightly change 
the definition of proton richness or neutron richness. As 
we can see from Fig. [21 superheavy nuclei have reached 
the /3-stability line around Z — 100 if we use the mi- 
croscopic /3-stability lines. However, there is still a gap 
if we use the phenomenological /3-stability line. In the 
following discussions, we will not distinguish the mi- 
croscopic /3-stability line and the phenomenological one. 
Second, all the three modern models predict quite similar 
most bound nuclei in each isotopic chain, which agree re- 
markably well with the experimental data. Third, the 
experimentally-synthesized superheavy nuclei populate 
to the left of the fission stability line and to the right 
of the /3-stability line, the same as what we observed for 
those heavy nuclei with 89 < Z < 100 H^. This finding 
leads us to speculate that the center of the superheavy is- 
land is most probably in between these two lines. We will 



call this region the "superheavy valley" . This speculation 
can find its support in "ordinary nuclei" , where we know 
most nuclei exist along the /3-stability line. But as we can 
see from Fig. |21 for them the fission stability line is quite 
close to the /3-stability line. For heavier nuclei, such as 
actinide nuclei, they populate in between (23. However, 
since we have not taken into account alpha-decay in this 
work, the exact location where the superheavy nucleus 
has the longest lifetime cannot be predicted. Lastly, the 
shell effects lead the most bound nucleus in each isotopic 
chain in the superheavy region to tend to have more pro- 
tons, compared to their light- or medium-mass counter- 
parts. The special role of neutron shell closure in this 
respect is most conspicuous at N — 50, 82 and 126. Af- 
ter three major neutron shell closures, the fission stability 
line deviates a lot from the /3-stability line in the super- 
heavy region, and thus, forms the so-called "superheavy 
valley" . 

The last observation above endows us with a powerful 
tool of identifying the next major neutron shell closure. 
That is to say there should also be a sudden increase in 
the fission stability line at the next major neutron shell 
closure. Using this argument, no neutron shell closure 
is found in the area we investigated, i.e. N < 180. It 
would be very interesting to study the N > 180 region, 
but it might be a difficult task due to the amount of com- 
puter resources needed. Fortunately, using our findings 
here, the searching region can be considerably reduced. 
Noticing that the fission stability line behaves as a linear 
function of N and Z between two major neutron shell 



5 



closures, we can approximate the theoretical and experi- 
mental fission stability lines (for N > 126) by the follow- 
ing formulae: 

Z = (17.105 ± 0.959) + (0.577 ± 0.006)iV (theory) (5) 
Z = 16.118 + 0.580A (Audi03). (6) 

The fission stability lines derived from the RMF+BCS, 
HFB-8 and FRDM calculations are all quite similar to 
each other, and therefore they have been averaged to ob- 
tain the above formula denoted by "theory" . One should 
note that these approximations are valid only up to the 
next major neutron shell closure. Future investigations, 
therefore, can be performed along these lines with the 
fission stability line as the upper bound. 

IV. SHAPES OF SUPERHEAVY NUCLEI 

Most experimentally-synthesized superheavy nuclei are 
believed to be deformed. This could be verified from two 
aspects: First, studies of "ordinary nuclei" revealed that 
only magic nuclei are spherical; second, actinide nuclei 
are known to be strongly deformed. In Fig. we plot the 
quadrupolc deformation parameters, /?2, of the 600 nuclei 
we calculated as functions of neutron number TV and pro- 
ton number Z. It is quite interesting to note that both 
RMF+BCS and HFB-8 calculations predict very strongly 
prolate deformations, > 0.45, for those nuclei in the 
upper-right corner. Even for those nuclei in the lower- 
left corner, the deformation is still appreciable, with (3% 
ranging from 0.15 to 0.35. However, there are two ma- 
jor differences between the predictions of the RMF+BCS 
calculations and those of the HFB-8 calculations. Firstly, 
there is a small number of spherical nuclei near N = 184 
and Z = 114 in the HFB-8 calculations, but there is no 
spherical nucleus at all in the RMF+BCS calculations. 
Second, strongly prolate deformation occurs in the HFB- 
8 calculations at mass numbers smaller than those pre- 
dicted by the RMF+BCS calculations; therefore it might 
indicate possible shape coexistence 0, HH, H3 • 

In other words, the RMF model with the effective force 
TMA predicts no doubly magic nucleus in the area we in- 
vestigated. The HFB-8 calculations seem to prefer iff 114 
to be a doubly magic nucleus, but the deformation pat- 
tern predicted for these superheavy nuclei are different 
from those observed for "ordinary nuclei" 0, , where 
a much larger number of nuclei near the doubly magic 
nucleus are found to be spherical. Therefore a decisive 
conclusion is not possible. This is, in fact, consistent 
with the recent study of Cwiok et al. j2jj, where they 
discussed the possibility of triaxial deformations for those 
nuclei near TV = 184 and Z = 114. Thus, this difference 
is once again the old model dependence problem. In Ref. 
, the spherical doubly magic nuclei are searched for all 
the usual parameterizations while here TMA is used to 
investigate their deformations. One should note that the 
octupole deformation and other beyond mean-field corre- 
lations can greatly reduce the second fission barrier |3lj| ; 
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FIG. 4: (color online). Proton quadrupole deformation pa- 
rameters, p2, of superheavy nuclei with Z = 101-120 and 
N — 134-190 as functions of neutron number N and proton 
number Z. The predictions of the RMF+BCS calculations 
are compared with those of the HFB-8 @| and FRDM [h} 
mass formulae. 

therefore, it would be interesting to investigate whether 
the conclusions of the present work would be modified if 
these correlations are taken into account. 

It is to be noted that the most remarkable difference 
exists between the microscopic models, RMF+BCS and 
HFB-8, and the macroscopic-microscopic model, FRDM: 
FRDM displays a transition from moderately deformed 
shapes to spherical shapes from the lower-left corner to 



6 



the upper-right corner. That is to say we find a similar 
pattern that we observed for those conventional doubly 
magic nuclei 0, |2£j, i.e. spherical shapes for doubly 
magic nuclei and those nearby. Therefore once again we 
conclude that f||ll4 is probably a doubly magic nucleus 
in the FRDM calculations, or more precisely, there is 
probably a major shell closure at N = 184. 

Of course, one should note that no clear indication of 
the N = 184 shell closure in the RMF+BCS and HFB- 
8 calculations (as shown in FigQJ docs not necessarily 
mean it does not exist. As well known nowadays, proton 
(neutron) shell closures are also neutron (proton) number 
dependent. Therefore, this implies if N = 184 is truly 
a neutron shell closure, we might have to look for the 
evidence in the region with Z > 120. Similarly, Z = 114 
or Z = 120 might be a major proton shell closure in the 
more neutron rich side. 



V. SUMMARY 

By performing a RMF+BCS calculation of about 600 
superheavy nuclei and employing the latest theoretical 



and experimental results, we have studied two extremely 
important subjects of superheavy nuclei, their stabilities 
and shapes. Both theory and experiment showed that all 
the experimentally-synthesized superheavy nuclei lie in 
between the fission stability line and the /3-stability line, 
i.e. the "superheavy valley" . It was also shown that the 
fission stability line and the /3-stability line tend to be 
more proton rich than their light- or medium-mass coun- 
terparts. In this sense, it is justified to say that the ob- 
served "proton richness" of superheavy nuclei is not only 
a result of the limitation of current experimental meth- 
ods but also a manifestation of their inherent nature. 
Although all the three theoretical models (RMF+BCS, 
HFB-8 and FRDM) predict most superheavy nuclei to 
be deformed, they differ from each other for those nuclei 
near N = 184. This model dependence might be removed 
by a more precise description of the single-particle spec- 
tra of the heaviest nuclei, which would be our next work. 
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